function quantify_comparsion()
load('weizmann_horse_db\sketch_result\horse_new.mat'); % our result
 results     =   gt_evaluat_result;
load('weizmann_horse_db\sketch_result\horse_comparison.mat'); % sketch-ucla
 results(:,:,2)     =   gt_evaluat_result;
 
load('dpm_horse_comparsion\result\horse_dpm_comparison.mat'); % pose-dpm
 results(:,:,3)     =   gt_evaluat_result;

load('HDT_horse\save\horse_HDT.mat'); % pose-dpm
is_valid = gt_evaluat_result>0;
gt_evaluat_result(~is_valid) = NaN;
 results(:,:,4)     =   gt_evaluat_result;
 
 is_valid(:,154)    =	[];
 results(:,154,:)   =   [];
 
 figure(1);  hold on;
 plot(results(1,:,1), 'r'); 
 plot(results(1,:,2), 'g');
 plot(results(1,:,3), 'b');
 plot(results(1,:,4), 'c');
 hold off;
 xlim([1 164]);
 xlabel('Sample No');  
 ylabel('Error');
set(gcf,'color','w');
legend('Our model','ABM', 'MPM', 'HDT');
% img  =  getframe;   imwrite(img.cdata, 'error.png');
 
[mean(results(1,:,1))  std(results(1,:,1))] 
[mean(results(1,:,2))  std(results(1,:,2))]
[mean(results(1,:,3))  std(results(1,:,3))]
[mean(results(1,is_valid(1,:),4))  std(results(1,is_valid(1,:),4))]
